Resonant excitations of single and two-qubit systems coupled to a tank circuit 
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The interaction of flux qubits with a low frequency tank circuit is studied. It is shown that 
changes in the state of the interacting qubits influence the effective inductance and resistance of the 
circuit, which is the essence of the so-called impedance measurement technique. The multiphoton 
resonant excitations in both single flux qubits and pairs of coupled flux qubits are investigated. In 
particular, we compare our theoretical results with recent spectroscopy measurements, Landau-Zener 
interferometry, and the multiphoton fringes. 
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I. INTRODUCTION 

Quantum properties of the atom-photon interaction 
have been mainly a subject of quantum optics and atomic 
physics. For example, the interaction between an atom 
and a photon was studied by experiments with a single 
atom in a cavity, this topic represents a textbook subject 
of quantum electrodynamics^ Recently a similar config- 
uration has been realized for superconducting quantum 
circuits^^ where Josephson qubits play the role of atoms, 
while a microwave resonator replaces the cavity. Usually 
for such kind of experiments the photon energy is close 
to the level separation A/h of the atom or qubit. In 
contrast, in the experiments of Refs. driven 
qubits are coupled to a LC tank circuit with a resonance 
frequency well below the level splitting of the qubits. In 
particular, our study here is motivated by the recent 
experiments^ where driven one or two flux qubits sys- 
tems were inductively coupled to a low frequency super- 
conducting tank circuit. 

The flux, or persistent current, qubit is a supercon- 
ducting ring with three Josephson junctions.— The cir- 
cuit is characterized by a two dimensional potential 
which, for suitable qubit and external parameters, ex- 
hibits two minima. If the applied magnetic flux is half 
a flux quantum, these minima have equal energies. In 
the flux-neighborhood of this point the degeneracy is 
lifted due to finite tunneling probability between min- 
ima. Therefore, the circuit forms an effective two-level 
quantum system. 

Strictly speaking a system of qubits coupled to a 
resonant circuit should be treated quantum mechan- 
ically, as in Refs. However, due 
to the weak coupling of the qubits to the classi- 
cal circuit the qubits-oscillator system can be treated 
semiclassicallyi^i^iii'i^ii^. In this paper we study the 
impact of the qubits on the tank circuit in terms of the 
effective inductance and resistance of the tank. We show 
that, in some limiting cases, the analyzed equations can 
be simplified resulting in a more transparent description 
of the behavior of the investigated qubit system. 



In this paper we will address the measurement with the 
resonant tank circuit in three different driving regimes of 
the qubits: in the ground state (without driving field), 
and in the weakly and strongly driven regimes. Such 
regimes are useful for controlling the state of the qubit 
with the driving field. Particularly interesting is the 
strongly driven regime, where due to the interference be- 
tween different Landau-Zener tunnelling events the state 
of the system quasi-periodically depends on both the 
DC bias and the AC driving amplitude3i2ii2^i2^iSi Moti- 
vated by our experimental paper—, we study the measure- 
ment technique in application to both single and coupled 
flux qubits. Systems of coupled qubits have been stud- 
ied previously^^ as well as spectroscopy with switching 
current readout i^i^l The study of the dynamical driv- 
ing of coupled qubits is important to form an entangled 
state (in particular, preparation of maximally entangled 
Bell states), and to perform two-qubit operations, such 
as a GNOT gatei^i^ Particularly, we study the multi- 
photon resonances in a two-qubit four-level system. This 
is analogous to the resonances studied in Ref. [13 where a 
multilevel system based on a single flux qubit including 
the upper levels was considered. 

The remainder of this paper is organized as follows. In 
Sec. II we derive equations which describe the influence 
of the qubits on the tank circuit in terms of the effective 
inductance and resistance. Limiting cases are considered 
in Sec. Ill for a single flux qubit coupled to the tank cir- 
cuit. In Sec. IV equations for the coupled qubits system 
are formulated. Numerically calculated results are pre- 
sented in Sees. V and VI for single and coupled qubits 
respectively. Experimental results for the multiphoton 
resonances in single flux qubit are shown in Sec. VII. 



II. DESCRIPTION OF THE MEASUREMENT 
WITH TANK CIRCUIT 

Consider a tank circuit which consists of an inductor 
Lt, a capacitor Ct, and a resistor Rt connected in par- 
allel (see Fig. [T]). The voltage V in the current-biased 
tank circuit (/bias = I A sinwrft), which is pierced by the 



2 



external flux , is described by the following non-linear 
equationJ^ 



where Q = ujtCtRt is the quality factor of the unloaded 
tank circuit (at $e = 0) and the functions (3s.c{v) are 
given by: 



CtV 



V 



V 



Rt L7 



(1) 



The external flux $e is assumed to be proportional to the 
coupling parameter and to depend on time via voltage 
V and its derivative V .^"^ 







■X — x- 





FIG. 1: Two-qubit system coupled to a tank circuit. The flux 
qubits are pierced by magnetic fluxes "l>i'''''' induced by the 
currents in the controlling coils (not shown in the scheme) 
and by the current in the tank's inductor. The qubits are 
coupled to each other and to the tank circuit. The resonant 
tank circuit consists of the inductor Lt, capacitor Ct, and 
resistor Rt; the circuit is biased with a RF current /bias. The 
tank voltage V is the measurable value. 

The stationary oscillations in the non-linear system de- 
scribed by Eq. ([T]) can be reduced to oscillations in a 
linear system by making use of the Krylov-Bogolyubov 
technique of asymptotic expansion. ■^^ Specifically, in the 
first order approximation with respect to the coupling 
parameter k'^ and close to the principal resonance (wrf « 
CUT = 1/\/LtCt) the equivalent linear system is char- 
acterized by the effective resistance i?off and inductance 
Leff as following-!^ 



27r 

PsM \ 1 /" i / , ■ ,^ f sini/' 1 , , 

o ) \ } ~ — / $e(wcosw, — wwrf smw) < , > . 
I3c{v) } nv J ^ ^' ^'{cosipj ^ 



(6) 

Let us also introduce notations >£'s,c for the partial deriva- 
tives: 

d<i>„ ■ (9$ ■■ 

= Xs{v,Tp)v simp + Kc{v,ip)''^ (^os^p, (7) 

where ip = uj^ft + a. If >Cs^c{v,ip) are functions of ei- 
ther si'mp or cos^p only, i.e. >Cs^ci'v,tp) = >^s c(''^i sin?/;) or 
>Cs^c{v,ip) = >c'^^^{v,cosip), the Eq. ([6]) leads to a simpler 
form for P^^c- 



>Cs{v, ip) SVC? ipdip, 



>Cc{v, ^p) cos^ Ipdip. 



(8) 



(9) 



As a result of Eqs. ^ and [5] the resonant frequency ojcs 
becomes amplitude-dependent and is shifted by: 



(10) 



The phase shift a and amplitude v depend on the fre- 



quency detuning ^0 



and the qubit state (via 



/3s. c). In the stationary regime they are given by: 



V = IaRcS cos a. 



tana = 2gfa/^, + i^j, 



L'T — L/Q 



(11) 



which can be rewritten in terms of the effective quality 
factor Qcff and effective frequency shift ^off : 



CtV 



^ + ^ 

RoS LcS 



V — V cos(ti>i-ft + a), 



Res{v) Rt 



Loff(w) Lt 



(2) 



(3) 



(4) 



(5) 



tana = 2(5eff^eff , 



Qos{v) = ujtCtRcs{v) 



iM^) = + ^fU^) = 



Q 



1-Qf3s{v)' 

LOT 



(12) 
(13) 
(14) 



Thus, the observable values - the amplitude v and the 
phase shift a - are defined by equations pT|) or (fT^ . 
which depend on the response of the measurable system. 
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III. MEASUREMENT OF THE PERSISTENT 
CURRENT QUBIT 

In this section we consider the system of the tank cir- 
cuit coupled to a persistent current qubit with geometri- 
cal inductance L and average persistent current /qb- The 
qubit is considered to be weakly coupled to the tank cir- 
cuit via a mutual inductance M . As we discussed above, 
strictly speaking, the dynamics of the tank circuit has 
to be considered jointly with the dynamics of the qubit. 
However, in this section we consider two limiting cases, 
when the dynamics of the qubit can be treated separately 
from the dynamics of the tank circuit. For simplification 
we introduce phenomenologically the relaxation time T\ 
(which can be caused by the tank as well) and consider 
the weak coupling limit, = -^j— <C 1. 



A. Low-quality qubit (Ti <C T): phase shift probes 
effective inductance of qubit 

When all the qubit 's characteristic times, and in par- 
ticular the relaxation time Ti , are smaller than the tank's 
period T = 2'k/ujt, the equations can be simplified, since 
the equations for the tank voltage can be averaged over 
the fast oscillating terms. This averaging is assumed 
henceforth for all values in this subsection. Then the 
time derivative of the flux $e, induced by the qubit in 
the tank circuit can be described as: 



MI, 



qb 



9$ 



(15) 



where $ = $dc + MIl is the flux in the qubit 's loop 
(herein we ignore the small self-induced fiux — L/qb), 
which consists of the time-independent part $dc and of 
the flux, induced by the current II in the tank's induc- 
tor. This can be rewritten by introducing the (inverse) 

-1 _ '9-fqb(t') 



effective inductance of the qubit, £ 



and the 



inductance value, which characterizes the response of the 
qubit, L = APC-\ Then it follows: 

$e = i(/L)/L; (16) 
and for the tank voltage we havei^ 

V^LTiL-^e^iLT~L{lL))iL- (17) 

Since L cx fc^ and L <g; Lt, in the first approximation in 
fc^ we can insert II in the r.h.s. of Eq. found from 
the above equation: 



hit) -Y^Jvdt 



sin(ijji.fi -I- a). 



lutLt 

Then by Eqs. fTlig)) we have 

j3, = 0, 13^ = 1- J k^LC-^{v, tP) cos^ T/'dV', 





(18) 



(19) 



where the qubit 's effective inductance is defined by the 
total flux <i>, piercing the qubit's loop: 



5$ 



-V sin ^ 



It follows: 



tana « 2QS,a + Ql3c 



IaRt cos a. 



(20) 



(21) 



Actually, this is a generalization of the result of Ref. [Tj 
for the case when the qubit can be in the excited state 
(which is taken into account by the expectation value of 
the current /qb)- 

If the bias current amplitude I a is small enough to be 
ignored in Eq. (I^Dl) . where v ^ IaRt, then:^^ 



tana 
£-1 



2Q^o + k^Q 

9/qb(^dc) 



L 
£' 



IaRt cos a, 



(22) 



which means that at the resonant frequency (^o = 0) the 
tank's phase shift a is approximately proportional to the 
inverse inductance of the qubit C~^. 



B. Higher-quality qubit (Ti < T): effective 
resistance due to qubit's lagging 

Consider the case when the qubit relaxation time Ti 
is of the same order as the tank's period T (namely, 
Ti < T). We assume exponential delayed response of 
the qubit 1 — exp(— t/Ti) which can be phenomenologi- 
cally described as a qubit response lagging relative to the 
tank circuit (as e.g. in Ref. ^3^. Therefore, Eq. is 
replaced by: 



$e(i) = i(/L(t')KL(0; 



(23) 



where t' — t — Ti is retarded time. Thus, the qubit's 
response depends on the current in the tank I^ = I Lit'), 
which is given by: 



lL{t') 



lutLt 

V 



sin(cjrft' + a) 

[C sin(u;rft + a) — S'cos(ajrf^ + en)) 



(24) 

where S — sin(ijJi.fTi) and C — cos(ijJifri). For the sake of 
simplicity we consider small bias current approximation, 
in this case Eqs. (|23II24I) and definitions result in 

gets: 



k^^S and /3c 



fc^^C. Then, from Eqs. (HO one 



Rt 



L 



—^Kl-S-k'Q-, 



R, 



£' 



(25) 
(26) 



4 



tan a 



2Q^o + C ■ eQL/C 

1-S-k^QL/C ' 
cos a 



IaRt l-S-k^QL/C 



(27) 
(28) 



Consider these expressions in the first approximation in 
k^QL/L] for = we obtain the following important 
result: 



tana ~ C ■ k^QL/C, 

1 + S-k^QL/C, 



(29) 



IaRt 



which, in particular, shows that as C — > there can 
be changes in the amplitude v without changes in the 
phase shift a as well as in shift of the resonant frequency 
(Eq. (fTU]) ). It is important to note that both the phase 
shift and amplitude are related to the qubit's effective 
inductance C, which explains their similar behavior in 
experiment. These equations might be useful for quali- 
tative analysis of experimental results. 



IV. COUPLED FLUX QUBITS 

A. Hamiltonian 

Close to its degeneracy point the flux qubit^° can be 
described by the pseudospin Hamiltonian: 



^^Iqb — 



-0-1 



eit) 



-CT3, 



(30) 



where the diagonal term e is the energy bias, the off- 
diagonal term A is the tunneling amplitude between the 
wells (which corresponds to the definite directions of the 
current in the loop) and aj are Pauli matrices. 

For the system of coupled qubits the effective Hamil- 
tonian is: 



2qbs 



E 

i=l,2 



J 



(31) 



where J is the coupling energy between qubits, and a{ , 
CTg*'' are the Pauli matrices in the basis , IT)} of the 
current operator in the i-th qubit; namely, a^P — aa'^do, 

(2) 

Ca — cro aa, <Jo is the unity matrix. The current op- 
erator is given by: li = —Ip''a^^\ with /p*^ the absolute 
value of the persistent current in the z-th qubit; then the 
eigenstates of correspond to the clockwise (0-3 \ l) = 
— IJ,)) and counterclockwise (0-3 \1) — current in the 
i-th qubit. The tunneling amplitudes are assumed 

to be constants. The biases ei — 2/p*''$o/*-*'' (^) ^^'^ con- 
trolled by the dimensionless magnetic fluxes /^*''(i) ~ 
^i/^o — 1/2 through i-th qubit. These fluxes consist of 
three components: 



/''nO=/. + ^ + /acSinu;t. 



(32) 



Here fi is the adiabatically changing magnetic flux, ex- 
perimentally applied by the coil and additional DC lines. 
The second term describes the flux induced by the cur- 
rent II in the tank coil, to which the z-th qubit is coupled 
with the mutual inductance Mi. And /acSino^t is the 
harmonic time-dependent component driving the qubit, 
typically applied by an on-chip microwave antenna. 



B. Entanglement 

It is convenient to present the density matrix for two 
qubits in the following form: 



Ral3 

Rqo 

Rob 



(To (8) CTo 4 

CTq (X) CTb 4 



RaO 
Rah 



(33) 



which was shown to be suitable for both the defini- 
tion and the calculation of the entanglement and other 
characteristics in multi-qubit system. ^^'^■^'■^^'^^ Here < 
{<^i /?} 5[ 3 and 1 < {a, b} < 3; the summation over 
twice repeating indices is assumed. The two vectors Rao 
and Rob, so-called coherence vectors or Bloch vectors, 
determine the properties of the individual qubits, while 
the tensor Rab (the correlation tensor) accounts for the 
correlationsj2S (In the notations of Ref. HI: Aa(l) = RaO, 
Xb{2) = Rob, Kab = Rab ) Following Ref. [H, we choose 
the measure of entanglement £ to be the following: 



3 ^ ' 



Rab — RaoR^ 



Qb- 



(34) 



This measure of entanglement fulfills certain reasonable 
requirements (0 < ^ < l)j in particular, £ — for any 
product state and £ = 1 for any pure state with vanishing 
Bloch vectors Rao and iJofci corresponding to maximum 
entangled states (see Ref. |3j for more detail) . 



C. Liouville equation 

The dynamics of the density matrix without taking 
into account the relaxation processes can be described 
by the Liouville equation: ip = [H,p\, which is gener- 
ally speaking a complex equation. We set both h ^ 1 
and ks = 1 throughout. By the proper choice of the 
parametrization the Liouville equation can be written in 
the form of a system with a minimal number of real equa- 
tions (3 - for one qubit and 15 for two qubits). 

To deal with the Liouville equation, we make use of 
the parametrization (decomposition) of the density ma- 
trix as described by Eq. (|33|) . This allows to benefit from 
the properties of the density matrix, namely from its her- 
miticity (then i?^^ are real numbers) and from the nor- 
malization condition, Trp = 1 (then i?oo = 1)- It follows 
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that the density matrix is parameterized by 15 indepen- 
dent real values. It is useful to note that: 



Tvipa^^A Rao, Tripa'^;' ) = Rot 



(35) 



After straightforward algebra the Liouville equation 
yields: 



RiO 


— F ■ R(^) 7? n - 


^ ^3ni J Rn3 i 




Roj 


— ^mnj R^yJ Ron 


+- £3nj J Rzn , 




R.,j 


- £ -B^^^R - 








+ Sjse^niJRnO 


- ^iS^SnjJRon- 


(36) 



where B^'^ are "local magnetic fields", which for the flux 
qubits are defined as B^'^ = (— Ai,0, — e^), and Smni is 
the Levi-Civita symbol. 



D. Effective inductance of coupled qubits 

For describing the effective inductance of coupled 
qubits it is important to note that the current in the 
z-th qubit depends on the fluxes in both qubits, /^j] = 

/qb(*i"\$i''^)- Here the fluxes consist of a DC 

part, $i, and of the flux generated by the current in the 

tank coil, 7^: ^i*' = + AIIl- For simplicity we as- 
sume here that the qubit-tank mutual inductance M is 
the same for both qubits (for more detail see Ref. [l9h . 
Then the time derivative is: 



f 9 



d 



(37) 



In the limit of small bias current in the tank we can 
substitute for <f>x and define the qubit effective in- 
ductance with the symmetric change of the flux bias in 
both qubits^e ^ ($^, ^^,)■. 



L7^ = 



d 



rii) 



( d 



7^)^^^- (38) 



Then for the case of low-quality qubits, when their char- 
acteristic times are smaller than the tank's period, analo- 
gously to Sec. III. A, we obtain at the resonance frequency 
(Co - 0): 



tana 



where we introduced the notation 



(39) 



(40) 



$0 

V. RESULTS FOR SINGLE FLUX QUBIT 

A. Spectroscopy 

In Ref. H it was shown that the qubit parameters can 
be determined both by measurements in the ground state 



or by employing a spectroscopic measurements when the 
qubit is resonantly excited. In this section we show re- 
lated numerically calculated graphs, making use of the 
results of previous sections. 

Consider a qubit biased with a DC flux <I>dc and driven 
with an AC flux <I>acSinct;t, introducing 

/dc = *dc/*0 - 1/2 and /ac = $ac/$0- 

In order to get the effective inductance £, as defined by 
Eq. (|20|) . we have to calculate the average current in 
qubit: /qb = (/) = Tr(p/), where / = —IpCr^ is the 
current operator defined with the amplitude Ip and the 
Pauli matrix a^. We calculate the reduced density ma- 
trix p with the Bloch equations, ''■^^'■^^ which include phe- 
nomcnological relaxation times, Ti and T2. It is conve- 
nient to express the density matrix in the energy repre- 
sentation: p = (1/2) (tq -t- Xti + Yt2 + Ztz), where Ti 
are the Pauli matrices for this basis and tq stands for the 
unity matrix. Z is equal to the difference between the 
populations of the ground and excited states. As a result 
the effective inductance is given byji^ 



d ( A 



2/„$ 



(41) 



where AE = y/A^ + (2/p$) . 

First consider the ground-state measurement, which is 
described by Eqs. (|19ll2ip . These equations not only al- 
low us to reproduce the results of Refs. [13 and [H, valid 
for the case where the system is in the ground state, 
but also describe the situation when the qubit is excited. 
Consider the influence of temperature, when the qubit is 
in a thermal mixture of the ground and excited states. 
In this case X = and Z = tanh(Ai;/2r). The result- 
ing tank phase shift is shown in Fig. [2] for the following 
parameters:?^ A/h = 2 • 0.65 GHz, Ip<i>o/h = 930 GHz, 
ujt/2tt = 32.675 MHz, L/p/$o = 0.0055, M/L = 0.725, 
Q = 725, k = 0.02. The accurate account of Z allows us 
to describe the widening of the phase shift dip, as shown 
in the inset in Fig. [21 which was reported in Ref. 39. The 
widening is due to the term that comes from differenti- 
ating the tanh in Eq. (j4ip : this term becomes relevant 
for temperatures larger than A = 1.3 GHz, and results 
in the exponential rise of the width for T > T* = A 
(^tanhx ~ 4exp(— 2a:) at x > 1). 

Now consider the spectroscopical measurement, where 
the qubit is driven with the AC flux. In Fig. we 
demonstrate the dependence of the phase shift a on the 
bias flux /dc at tOrf — l^jt for different driving amplitudes 
at the driving frequency uj/2Tr = 4.15 GHz. The param- 
eters for plotting the graph were taken as in the related 
experiment^ and are given in Table [D In Fig. [3] the up- 
per curves are shifted vertically for clarity. In Fig. ^h) 
we plot the amplitude v versus the bias flux /dc with 
the phenomenological lagging parameter S = 0.8 for sev- 
eral values of the driving frequency w, making use of 
Eqs. (|27ll28p . In the experimental case the positions of 
these resonances at a given driving frequency allow to 
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FIG. 2: (Color online). Influence of temperature on the 
ground state measurement: the dependence of the tank phase 
shift on the flux detuning fdc = $dc/$o — 1/2, when the qubit 
is thermally excited. The numbers next to the curves stand 
for the temperature in GHz. Inset: temperature dependence 
of the width A/dc of the dip at half-depth in the phase shift, 
shown in the main panel. 

TABLE I: The parameters used for plotting Fig. O A, /p<&o, 
Fi, and r2 are tunneling amplitude, energy bias, relaxation 
rate, and dephasing rate, respectively. The H describes the 
coupling between qubit and tank circuit. All parameters are 
in units of /i-GHz except H which is dimensionless. 



A 




r 


Fi 


Fa 




3.5 


700 


1.4 


0.7 


0.7 


2.6 • 10"^ 



determine the energy structure of the measured qubit 
Figure [3] demonstrates the effect described in section III: 
for S both the phase shift a and the ampUtude v 
depend on C~^, which results in the alternation of peak 
and dip around the location of the resonances. We also 
note that, in addition, the terms nonlinear in S distort 
this structure. The second-order term for example, is 
proportional to £"^(35^ — 1). For values S > l/V^ this 
term leads to increasing of the peak and decreasing of the 
dip leading to the asymmetry visible in the upper curve 
of Fig.^b). 

In Fig. [4] we plot the phase shift a and the amplitude 
V as functions of the bias current frequency uj^f and the 
flux detuning fdc with the phenomenological lagging pa- 
rameter S for the flux qubit with the parameters given in 
Table. [J The dashed white line shows the tank resonance 
frequency uj^{/2tt — lut/'^t^ = 20.8 MHz. Note that for 
the lagging parameter close to 1 (5* = 0.8) the changes 
in the phase shift in Fig. SJa) are small at the resonance 
frequency (along the white line) while the voltage am- 
plitude in Fig. mjb) changes substantially. Such changes 
of the tank effective resistance or, equivalently, quality 
factor were studied in Ref. [^. 




FIG. 3: (Color online). Curves for a resonantly excited flux 
qubit: (a) the phase shift and (b) the amplitude of the tank 
voltage versus flux detuning fdc- In panel (a) the numbers 
next to the curves stand for the driving amplitude /ac multi- 
plied by 10^; in panel (b) the curves are for the amplitudes 
/ac • 10'^ — 1, 1.5, 3 from bottom to top. 



B. Landau-Zener interferometry 

LZ interferometry is demonstrated in Fig. [5] as the de- 
pendence of the tank voltage phase shift a on the mi- 
crowave amplitude /ac and the DC flux bias fdc- The 
qubit parameters were taken as for Fig. [3] (Table |T| and 
uj/2n = 4.15 GHz. 

The multiphoton resonances at discrete DC bias fdc 
(which controls the distance between energy levels) are 
clearly visible. These resonances appear when the energy 
of n photons matches the qubit's energy levels: n ■ fku = 
AE{fdc)- Note also the quasi-periodical character of the 
dependence on the AC flux amplitude /ac , this is known 
as Stiickelberg oscillations; the comparison of such graph 
to the experimental analogue^ (namely the estimation of 
the period of Stiickelberg oscillations) allows the relation 
of the microwave power to the AC flux amplitude /ac to 
be determined. 
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FIG. 4: (Color online). Multi-photon excitations of a single 
flux qubit. (a) dependence of the phase shift a and (b) the 
amplitude v on the bias current frequency uj^f and the flux 
detuning /dc; parameters: /ac = 8 ■ 10~^, uj/2-k = 4.15 GHz, 
S = 0.8. 



FIG. 5: (Color online). Landau-Zener interferometry: de- 
pendence of the tank phase shift on the flux detuning /dc = 
^dc/^o — 1/2 and on the driving flux amplitude /ac = '^ac/^o- 




0.02 



C. Impact of finite bias current 

Consider the impact of the finite bias current on the 
tank's response. When the bias current is small, its influ- 
ence can be neglected and there is a peak-and-dip struc- 
ture, around the point where the qubit is resonantly ex- 
cited to the upper state, see Eq. (|22p : these structures 
are distorted by the non-linear terms when the current is 
increased. This is demonstrated in Fig. [SI where we plot 
the dependence of the phase shift a on the bias flux /dc 
for two different values of the bias current amplitude Ia- 
The qubit parameters were taken as for Fig. [3] (Table. 



FIG. 6: (Color online). Impact of the flnite bias current: de- 
pendence of the phase shift a on the bias flux /dc for different 
values of the bias current amplitude I a, for uj/2-n = 4 GHz, 
/ac = 8 ■ 10-^ 



in terms of its energy structure, entanglement measure, 
and the observable tank circuit phase shift. Then we will 
study the one-photon excitation of the system, which is 
used for the spectroscopic measurements. Finally we will 
demonstrate the multiphoton fringes. For calculations we 
will use the two sets of parameters as in Refs.[25l (i) andU 
(ii) which are given in Table [Til Here, for clarity, we con- 
sider the case when the characteristic measurement time 



VI. RESULTS FOR THE TWO-QUBIT SYSTEM 

A. Resonant excitation 

In this section we investigate the resonant excitation of 
a system of two coupled flux qubits. First we will describe 
the effects of resonant excitation in a system of two qubits 



TABLE II: The parameters of the two-qubit systems. All 
parameters are in units of ^-GHz except H^.b which are di- 
mensionless. 





A, 


Ai, 






J 




' — 'a 


—6 


(i) 


1.1 


0.9 


990 


990 


0.84 


1.8 


■ 10"=* 


1.8 • 10"^ 


(ii) 


15.8 


3.5 


375 


700 


3.80 


1.4 


■ 10"=^ 


2.6 • 10"^ 
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T is larger than the characteristic times of the dynamics 
of the qubit (Ti ) . Then the tank circuit actually probes 
the incoherent mixture of qubit 's states and the time- 
averaged values of phase shift and entanglement should 
be considered. We calculate the energy levels (by diago- 
nalizing the stationary Hamiltonian) , the density matrix 
p, the observable tank circuit phase shift a (which is de- 
fined with the effective inductance of the qubits), and 
the entanglement measure E by making use of equations 
P4j|39p . Then we plot Fig. [7] for the set of parameters (i) 
and the driving frequency Lojl-K = 4 GHz, assuming the 
symmetrical change of the DC flux: fa — fb = fdc- Four 
energy levels are plotted in Fig. [Tlja). When the energy 
of n photons {n ■ huj) matches the energy difference be- 
tween any two levels Ej and Ei, the resonant excitation 
to the upper level is expected. Respectively, with the 



green (gray), black, and magenta (dark gray) arrows of 
the length 4, 8, and 12 GHz we show the places of possi- 
ble one-, two-, and three-photon excitations. The time- 
averaged total probability of the currents in two qubits 
to flow clockwise, Z — i?03 + is shown in Fig. [Tjb) 
to experience resonant excitation; black and red (dark 
gray) hues correspond to /ac — and 10~^. The reso- 
nances appear as hyperbolic-like structures in the phase 
shift dependence in Fig. [Tljc). The time-averaged entan- 
glement measure £ has a peak in resonance, Fig. El^d). 
The entanglement measure in a resonance increases due 
to the resonant formation of the superposition of states; 
this provides a tool for resonantly controlling the entan- 
glement and on the other hand comparing Figs. [71(c) and 
[7fd) is a method to probe the entanglement. 




FIG. 7: (Color online). Multiphoton excitations of coupled flux qubits: (a) energy levels, (b) total probability of the currents 
in qubits to flow clockwise Z, (c) the tank phase shift a, (d) the entanglement measure £ ~ versus the bias fdc ~ fa = ft- 



B. Spectroscopy (one-photon excitation) 

A weakly driven system of qubits can be resonantly 
excited to an upper level when the driving frequency 
matches the distance between the energy levels. This 
can be used for the spectroscopic measurement of the 
energy structure of the system. In analogy to the re- 
lated experiment^ we plot in Fig. [8] the tank phase shift, 
when the tank is coupled to the two flux qubits with the 
set of parameters (ii) and for the driving frequency and 
amplitude: uj/2n = 17.625 GHz and /ac = 10"^. The 



resonances are visualized as the ridge-trough lines, as de- 
scribed above. The wide dip around /& ~ is due to 
the ground state curvature. Solid lines are the energy 
contour lines which show where the photon energy hu 
matches the respective levels: ground and first excited 
(black), ground and second excited (white), and second 
excited and upper (third excited) levels (orange/gray). 

An interesting situation arises when one photon ex- 
cites the system to an intermediate level (to the second 
excited level in the figure), then, due to the non-zero 
level population, another photon can excite the system 



to a higher level (to the third excited level in the fig- 
ure). This is shown with the white circle in Fig. [51 which 
marks the region where the signal is increased. The grey 
circle shows that such increase may be not visible, when 
the former excitation (tn the second excited level") cre- 
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FIG. 8: (Color online). Resonant excitation of two coupled 
flux qubits: the tank phase shift versus partial bias fluxes in 
two qubits, fa and /;,. 



C. Multiphoton excitations 
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FIG. 9; (Color online). Multiphoton resonances in coupled 
qubits: dependence of the tank phase shift on the flux detun- 
ing in qubits a and 6, fa and fb- 



When the qubits are strongly driven, excitations due to 
multiphoton processes become possible. In this case reso- 
nances appear where the energy level difference is a mul- 
tiple of the photon energy fku. The multiphoton fringes 
are shown in Fig. [5] for the following driving and qubit 
parameters: (a) lu/2tt — 4 GHz, /ac = 5 • 10^'^, and qubit 
parameters (i); (b) uj /2t: = 4.15 GHz and /^c = 7 • 10~^, 
and qubit parameters from (ii) with equal parameters 
(those for qubit a are taken to be the same as for qubit 
b: Aa = Ab = 3.5, etc.). 

The lines in Fig. [9^b) are the energy contour lines. 
They show that the multiphoton resonances are mostly 
due to the excitation to the first excited level (black lines) 
with the interruptions in ridge-trough resonant structures 
(change of the signal), where higher levels are matched 
with the multiple photon energy; possible excitations to 
the second excited level are shown with the magenta 
(dash gray) line and to the upper level with the red (solid 
dark gray) line. 

VII. RELATION TO EXPERIMENT 

As mentioned above, experimental results related to 
Figs. 2, 3b, and 5 were already reported in Refs. [sol 



and 8. The presented theory gives good agreement with 
those experiments. Here we present some additional data 
obtained by making use of the same experimental proce- 
dure as in Ref. S In these experiments the tank circuit 
was realized as a parallel connection of a superconduct- 
ing niobium coil and commercial ceramic capacitor (see 
Fig. 1). The superconducting persistent current qubits, 
three Josephson junctions closed into a superconducting 
ring,^° were deposited by an aluminum shadow evapora- 
tion technique in the center of the superconducting coil. 
The voltage across the tank circuit V = vexp(ia) was 
amplified by a cryogenic amplifier and measured by an 
rf lock-in amplifier. 

In Fig. [To] multi-photon resonances for a single qubit 
are shown in dependence of both the DC flux f^c and the 
frequency of the driving current Wrf • This data should be 
compared with the theoretical result presented in Fig [Jj 
Additional multi-photon lines, similar to the theoretical 
results of Fig. [3^, are shown in Fig. [TTb . Those results 
were recorded at the resonant frequency only. Finally we 
studied the impact of finite bias current on the tank's re- 
sponse as shown in Fig. Illb which should be compared 
with Fig. [S] Also in this more detailed comparison with 
the experiments we find a good agreement with the the- 
ory presented in this work. Experimental results for the 
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multi-photon resonances in coupled qubits will be pub- 
lished elsewhere. 
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FIG. 10: (Color online). Experimental demonstration of 
multi-photon excitations in a single flux qubit. (a) Depen- 
dence of the phase shift a and (b) the amplitude v on the bias 
current frequency oj^f and the flux detuning fdc- Recorded for 
a microwave excitation with uj/2n = 4.15 GHz and amplitude 
/ac =4.5 ■ 10-^ 



FIG. 11: (Color online). Experimental dependence of the 
phase shift a on the bias flux fdc'- (a) multiphoton reso- 
nances for driving amplitudes /ac/10~"^ — 0.5, 2.5, 4.5, and 
6.5 from top to bottom; (b) influence of the bias current (red 



line /aM/c&o = 2-10" 
for lo/2tt = 4.15 GHz. 



and blue line - /aM/$o = 8 ■ 10"^) 



VIII. CONCLUSION 

The application of the impedance measurement tech- 
nique for monitoring resonant excitations in flux qubits 
was studied. It was shown that excitation of the qubits 
change their effective inductance. This results in a 
changed effective inductance and resistance of the tank 
circuit. Thus, the observable quantities - the amplitude 
and phase shift of the tank voltage - reffect the changes 
of the effective inductances of the qubits. Transparent 
expressions were derived in two limiting cases, one where 
the dynamics of the tank is slow relatively to the time 
scale of the qubits and another where both timescales 
are of the same order. In the first case the changes in 
the tank's resistance are negligible and the phase shift is 
directly related to the effective inductance. It has also 
been demonstrated that in the latter case both the ef- 
fective inductance and resistance of the tank are defined 



by the inductance of the qubit. Our theoretical analysis 
describes the state of the system of qubits in terms of the 
tank's effective inductance and resistance. This allowed 
us to describe the experimental results on multi-photon 
resonant excitation of single and coupled qubits. 

The impact of the finite bias current and of the temper- 
ature was described as well. In particular it was shown 
that the thermal excitation of the qubit to the upper level 
results in the widening of the dip visible in the phase shift 
as was observed previously in experiments^. 

The resonant excitation of qubits can be used not only 
for controlling their state, but also to determine their 
parameters. The developed theory has been used to re- 
produce the experimental results presented in Ref. H for 
single qubits and pairs of coupled fiux qubits. In ad- 
dition to single-photon ("spectroscopic") resonances we 
have shown the appearance of the multi-photon excita- 
tions in flux qubits coupled to the tank circuit. 
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